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Abstract 

A review of different cosmological models in diverse dimensions leading to a relatively small time 
variation of the effective gravitational constant G is presented. Among them: 4-dimensional general 
scalar-tensor model, multidimensional vacuum model with two curved Einstein spaces, multidimensional 
model with multicomponent anisotropic "perfect fluid", S-brane model with scalar fields and two form 
field etc. ft is shown, that there exist different possible ways of explanation of relatively small time 
variation of the effective gravitational constant G compatible with present cosmological data (e.g. ac- 
celeration): 4-dimensional scalar-tensor theories or multidimensional cosmological models with different 
matter sources. The experimental bounds on G may be satisfied ether in some restricted interval or for 
all allowed values of the synchronous time variable. 
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1 Introduction 



In the development of rclativistic gravitation and dynamical cosmology after A. Einstein and A. Friedmann 
there are three distinct stages: first, investigation of models with matter sources in the form of a perfect fluid, 
as was originally done by Einstein and Friedmann. Second, studies of models with sources as different physical 
fields, starting from electromagnetic and scalar ones, both in classical and quantum cases (see [T]). And 
third, which is really topical now, application of ideas and results of unified models for treating fundamental 
problems of cosmology and black hole physics, especially in high energy regimes and for explanation of the 
greatest challenge to modern physics - the present acceleration of the Universe, dark matter and dark energy 
problems. Multidimensional gravitational models play an essential role in the latter approach. 

The necessity of studying multidimensional models of gravitation and cosmology [21 El [4] is motivated by 
several reasons. First, the main trend of modern physics is the unification of all known fundamental physical 
interactions: electromagnetic, weak, strong and gravitational ones. During the recent decades there has been 
a significant progress in unifying weak and electromagnetic interactions, some more modest achievements in 
GUT [5], super-symmetric, string and super-string theories. 

Now, theories with membranes, p-branes and more vague M-theory are being created and studied. Having 
no definite successful theory of unification now, it is desirable to study common features of these theories 
and their applications to solving basic problems of gravitation and cosmology. If we really believe in unified 
theories, the early stages of the Universe evolution and black hole physics, as unique superhigh energy 
regions and possibly even low energy stage, where we observe the present acceleration, are the most proper 
and natural arena for them. 

Second, multidimensional gravitational models, as well as scalar-tensor theories of gravity, are theoretical 
frameworks for describing possible temporal and range variations of fundamental physical constants [U [HI El 
[§]. These ideas have originated from papers of E. Milne (1935) and P. Dirac (1937) on relations between the 
phenomena of micro- and macro-worlds, and up till now they are under thorough study both theoretically 
and experimentally. The possible discoveries of the fine structure constant and proton to electron ratio 
variations are now at a critical further investigation. 

Lastly, applying multidimensional gravitational models to basic problems of modern cosmology and black 
hole physics, we hope to find answers to such long-standing problems as singular or nonsingular initial states, 
creation of the Universe, creation of matter and its entropy, cosmological constant, coincidence problem, 
origin of inflation and specific scalar fields which may be necessary for its realization, isotropization and 
graceful exit problems, stability and nature of fundamental constants [6l [9l [10l Qj] , possible number of extra 
dimensions, their stable compactification, new revolutionary data on present acceleration of the Universe, 
dark matter and dark energy [10] etc. 

Bearing in mind that multidimensional gravitational models are certain generalizations of general rela- 
tivity which is tested reliably for weak fields up to 0.0001 and partially in strong fields (binary pulsars), it is 
quite natural to inquire about their possible observational or experimental windows. From what we already 
know, among these windows are [91 110] : 

- possible deviations from the Newton and Coulomb laws, or new interactions, 

- possible variations of the effective gravitational constant with a time rate smaller than the Hubble one, 

- possible existence of monopole modes in gravitational waves, 

- different behavior of strong field objects, such as multidimensional black holes, wormholes and p-branes, 

- standard cosmological tests, 

- possible non-conservation of energy in strong field objects and accelerators, if brane- world ideas about 
gravity in the bulk turn out to be true etc. 

Since modern cosmology has already become a unique laboratory for testing standard unified models of 
physical interactions at energies that are far beyond the level of existing and future man-made accelerators 
and other installations on the Earth, there exists a possibility of using cosmological and astrophysical data 
for discriminating between future unified schemes. Data on possible time variations or possible deviations 
from the Newton law as a new important test should also contribute to the unified theory choice and viable 
cosmological model choice as well [TU [13] . 

As no accepted unified model exists, in our approach [21 HI [14] we adopted simple (but general from 
the point of view of number of dimensions) models, based on multidimensional Einstein equations with or 



2 



without sources of different nature: 

- cosmological constant, 

- perfect and viscous fluids, 

- scalar and electromagnetic fields, 

- their possible interactions, 

- dilaton and moduli fields with or without potentials, 

- fields of antisymmetric forms (related to p-branes) etc. 

Our program's main objective was and is to obtain exact self-consistent solutions (integrable models) 
for these models and then to analyze them in cosmological, spherically and axially symmetric cases. In our 
view this is a natural and most reliable way to study highly nonlinear systems. It is done mainly within 
Riemannian geometry. Some simple models in integrable Weyl geometry and with torsion were studied 
as well. In many cases we tried to single out models, which do not contradict available experimental or 
observational data on variations of G. In some cases we used our methods for arbitrary dimensions in 
studying 4D models also. 

As our model [2J [31 3] we use n Einstein spaces of constant curvature with sources as (m+l)-component 
perfect fluid, (or fields or form-fields,), cosmological or spherically symmetric metric, manifold as a direct 
product of factor-spaces of arbitrary dimensions. Then, in harmonic time guage we show that Einstein 
multidimensional equations are equivalent to Lagrange equations with non-diagonal in general mini-super- 
space metric and some exponential potential. After diagonalization of this metric we perform reduction to 
sigma-model [TSJ [TB] and Toda-like systems [T7], further to Liouville, Abel, generalized Emden- Fowler Eqs. 
etc. and try to find exact solutions. We suppose that behavior of extra spaces is the following: they are 
constant, or dynamically compactified, or like torus, or large, but with barriers, walls etc. 

So, we realized and continue to realize the program in arbitrary dimensions (from 1988) [21 [3j (H [TH [10] : 
In cosmology we obtained exact general solutions of multidimensional Einstein equations with sources: 

- A , A + scalar field (e.g. nonsingular, dynamically compactified, inflationary) [18] ; 

- perfect fluid, PF + scalar field (e.g. nonsingular, inflationary solutions) [H)l [2"U] : 

- viscous fluid (e.g. nonsingular, generation of mass and entropy, quintessence and coincidence in 2- 
component model); 

- stochastic behavior near the singularity, billiards in Lobachevsky space, D=ll is critical, ip destroys 
billiards (1994). 

(For all above cases Ricci-flat solutions were obtained for any n) 

- also solutions with curvature in one factor-space; 

- with curvatures in 2 factor-spaces only for total N=10, 11; 

- with fields: scalar, dilatons, forms of arbitrary rank [21] (1998) - inflationary, A generation by forms 
(p-branes) [22] : 

- first billiards for dilaton-forms (p-branes) interaction (1999); 

- quantum variants (solutions of WDW-equation [18 ] 120 ] 123 ] ) for all above cases where classical solutions 
were obtained; 

- dilatonic fields with potentials, billiard behavior for them. 

For many of these integrable models we calculated also the variation with time of the effective gravitational 
constant. Comparison with present experimental bounds allowed to choose particular models or single out 
some classes of solutions. 

Similar methods were used for obtaining exact solutions in spherical symmetry case [2J [3J [H [TH [TU] • 
Dirac's Large Numbers Hypothesis (LNH) is the origin of many theoretical studies of time-varying G. 
According to LNH, the value oiG/G should be approximately the Hubble rate. Although it has become clear 
in recent decades that the Hubble rate is too high to be compatible with experiments, the enduring legacy 
of Dirac's bold stroke is the acceptance by modern theories of non-zero values of G/G as being potentially 
consistent with physical reality. 

There are three problems related to G, whose origin lies mainly in predictions of unified models of physical 
interactions: 

1) absolute G measurements, 2) possible time variations of G, 3) possible range variations of G - non- 
Newtonian, or new interactions. 
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After the original Dirac hypothesis some new ones appeared and also some generalized theories of grav- 
itation admitting variations of the effective gravitational coupling. We can single out three stages in the 
development of this field [6] : 

1. Study of phenomenological theories and hypotheses with variations of FPC, their predictions and 
confrontation with experiments (1937-1977). 

2. Creation of theories admitting variations of an effective gravitational constant in a particular system 
of units, analyses of experimental and observational data within these theories pQ (1977-present). 

3. Analysis of FPC variations within unified models [3] (present). 

Different theoretical schemes lead to temporal variations of the effective gravitational constant [53] : 

- Empirical models and theories of Dirac type, where G is replaced by G(t). 

- Numerous scalar-tensor theories of Jordan-Brans-Dicke type where G depending on the scalar field <fi(t) 
appears. 

- Gravitational theories with a conformal scalar field arising in different approaches [25] [I] (they can be 
treated as special cases of scalar-tensor theories). 

- Multidimensional unified theories in which there are dilaton fields and effective scalar fields appearing 
in our 4-dimcnsional spacctime from additional dimensions [3] . They may help also in solving the problem of 
the variable cosmological constant (from Planckian to present values) and the cosmic coincidence problem. 

A striking feature of the present status of theoretical physics is that there is no satisfactory theory unifying 
all four known interactions; most modern unification theories do not admit unique and universal constant 
values of physical constants and of the Newtonian gravitational coupling constant G in particular. Although 
the bounds on G and G(r) are in some classes of theories rather wide on purely theoretical grounds since any 
theoretical model contains a number of adjustable parameters, we note that observational data concerning 
other phenomena, in particular cosmological data, may place limits on possible ranges of these adjustable 
parameters. But, in any case variations of G may be an additional test of unified models, generalized theories 
of gravitation and cosmological models as well [13] . 

Here we restrict ourselves to the problem of G (for G(r) see [3] [H] [T] [7] ) • We show that various theories 
predict the value of G/G to be 10~ 12 /yr or less. The significance of this fact for experimental and obser- 
vational determinations of the value of or upper bound on G is the following: any determination with error 
bounds significantly below 10~ 12 /yr (combined with experimental bounds on other parameters) will typically 
be compatible with only a small portion of existing theoretical models and will therefore cast serious doubts 
on the viability of all other models. In short, a tight bound on G, in conjunction with other astrophysical 
observations, will be a very effective "theory killer" and/or significantly reduce the class of viable theories. 
Any step forward in this direction will be of utmost significance [26] [27] [28] . 

Some estimations for G were done long ago in the frames of general scalar tensor theories using the values 
of cosmological parameters (tt, H, q etc) known at that time [23 EOl [3TJ [32l [3J [TJ . It is easy to show that for 
modern values they predict G/G at the level of 10 _13 /yr and less (see also estimations of A. Miyazaki 33J, 
predicting time variations of G at the level of 10 _13 yr _1 for the Machian-type cosmological solution in the 
Brans-Dickc theory, Y. Fujji (sec in 34 ), of J.P.Mbelek and M.Lanchieze-Ray [35] on the level of 10~ 17 /yr 
for the simple 5D KK-theory with an external scalar field and section 2 of the present paper). 

The most reliable experimental bounds on G/G (radar ranging of spacecraft and planets dynamics [36(137] 
and laser lunar ranging [38] [39]) give the limit less than 10~ 12 /yr, so any result at less than this level will 
be very important for solving the fundamental problem of variations of constants and for discriminating 
between viable unified theories and cosmological models. So, further data on lunar laser and radar ranging 
and realization of MICROSCOPE, ASTROD, LATOR space projects and such multipurpose new generation 
type space experiment as Satellite Energy Exchange (SEE) for measuring G, absolute value of G and Yukawa 
type forces at meters and Earth radius ranges [40] become extremely topical. 

In what follows, we shall discuss predictions for G from generalized scalar-tensor theories and some 
multidimensional models. 

2 Scalar- Tensor Cosmology and Variations of G 

The purpose of this section is to estimate the order of magnitude of the gravitational constant G variations 
due to cosmological expansion in the framework of general scalar-tensor theories (STT) of gravity [41] , 
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Consider the general (Bermann-Wagoner-Nordtvedt) class of STT where gravity is characterized by a 
metric g^ and a scalar field tfi; the action is 

S = J d 4 x^g{f(4,)TZ[g] + h^g^^^ - 2U(cj>) + L m }. (2.1) 

Here lZ[g] is the scalar curvature, g — | det(g /j;y )|; /, h and U are certain functions of (f>, varying from theory 
to theory, L m is a matter Lagrangian. 

This formulation of the theory corresponds to the Jordan conformal frame, in which matter particles 
move along geodesies and hence the weak equivalence principle is valid and non-gravitational fundamental 
constants do not change. In other words, this is the frame well describing the existing laboratory, geophysical 
and cosmological observations. 

Among the three functions of 4> entering into (|2.ip only two are independent since there is a freedom 
of transformations <fi = 0(0 n ew) [29]. We use this arbitrariness, choosing h(<fi) = 1, as is done, e.g., in [42] . 
Another standard parametrization is to put /(</>) = </> and h(4>) = u>(4>)/(j) (the Brans-Dicke parametrization 
of the general theory (|2.ip ) In our parametrization h=l, the B-D parameter lo{4>) = f{f<j>)~ 2 \ the subscript 
(f> denotes a derivative with respect to <fi. The B-D STT is the particular case ui = const, so that in (|2.ip 



/(</>) =0 2 /(M, h=l. (2.2) 

For the conformal scalar field case see [1] [29] . 
The field equations that follow from (|2.ip read 

U4>- ±Kf <t> + U< t , = 0, (2.3) 

f(4>)fe - \8pi) = -4>^ v + h^' a <f>,a - spite) + (v M v- - s;a)f- r; (m) , (2.4) 

where □ is the D'Alembert operator, and the last term in (|2.4p is the energy-momentum tensor of matter. 
Consider now isotropic cosmological models with the standard FRW metric 



ds 2 = dt 2 -a 2 (t) 



1 — kr A 



(2.5) 



where a(t) is the scale factor of the Universe, and k = 1, 0, —1 for closed, spatially flat and hyperbolic 
models, respectively. Accordingly, we assume <p = <p(t) and the energy-momentum tensor of matter in the 
perfect fluid form T^( m ) = diag(,o, — p, —p, —p) {p is a density and p is a pressure). 
The field equations in this case can be written as follows: 

+ 3-0- ^r{ad + d 2 + k) + U tj} = 0, (2.6) 
a or 

^a 2 + k)^^ 2 + U~3-f + p, (2.7) 

4(2aa + a 2 + fc) = -U 2 + U - f - 2-f - p. (2.8) 

To connect these equations with observations, let us fix the time t at the present epoch (i.e., consider 
the instantaneous values of all quantities) and introduce the standard observables: H = a/a (the Hubble 
parameter), q — —ad/ a 2 (the deceleration parameter), f2 m = p/ p cr (the matter density parameter), where 
p cr is the critical density, or, in our model, the r.h.s. of equation (|2.6p in case k — 0: p cr = 3fH 2 . This is 
slightly different from the usual definition p CI = 3H 2 /8irG. The point is that the locally measured Newtonian 
constant in STT differs from l/(87r/); provided the derivatives and are sufficiently small, one has 

m 

Since, according to the Solar-system experiments, u> > 40000, for our order-of-magnitude reasoning we can 
safely put 8irG = 1//, and, in particular, our definition of p cr now coincides with the standard one. 
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The time variation of G, to a good approximation, is 

G/G « -./// = gH, (2.10) 

where, for convenience, we have introduced the coefficient g expressing G/G in terms of H. 

Equations (12.6|) - (12.8|) contain too many arbitrary parameters for making a good estimate of g. Let us 
now introduce some restrictions according to the current state of observational cosmology: 

(i) k — (a spatially flat cosmological model, so that the total density of matter equals p CI ); 

(ii) p = (the pressure of ordinary matter is negligible compared to the energy density) ; 

(iii) p = 0.3 ,o cr (the ordinary matter, including its dark component, contributes to only 0.3 of the critical 
density; unusual matter, which is here represented by the scalar field, comprises the remaining 70 per cent). 

Then equations (|2.6[) and (|2.8[) can be rewritten in the form 

^ 2 + U-3Hf = 2.lH 2 f, -^ 2 + U-2Hf-f = (l-2q)H 2 f. (2.11) 

Subtracting second from the first one, we exclude the "cosmological constant" U and obtain 

<j> 2 -Hf + f = (l.l + 2q)H 2 f. (2.12) 
The first term in equation (I2.12p can be represented in the form 

4> 2 = f 2 {df/d4>r 2 = /V/, (2-i3) 

and /// can be replaced with —gH . The term / can be neglected for our estimation purposes for an arbitrary 
function / and potential U(<j>). 

Then, (I2.12| divided by H 2 f leads to the quadratic equation with respect to g: 

ug 2 + g- q' = 0,q' = 1.1 + 2q. (2.14) 

According to present observations, the Universe is expanding with an acceleration, so that the parameter q 
is, roughly, —0.5 ± 0.2, hence we can take |g'| < 0.4. 

In case q' = we simply obtain g = —1/lu. Assuming 

H = hwo ■ 100 km/(s.Mpc) » h 1Q0 ■ 10~ 10 yr" 1 (2.15) 

and present limit u> > 40000, we come to the estimate 

\G/G\ < 4-10- 15 / lloo yr~\ (2.16) 

where h\oo is, by modern views, close to 0.7. So (|2.16|) becomes 

\G/G\ < 4-10- 15 yr- 1 . (2.17) 

For nonzero values of q' , solving the quadratic equation (|2 . 14[) and assuming q'uo 3> 1, we arrive at the 
estimate \g\ ~ y/q' /uj, so that, taking q 1 = 0.4 and again u > 40000, we have instead of (|2.16|) 

\G/G\ < 0.9-10- 13 /iioo yr _1 ~ 0.7-10 _13 yr _X , (2-18) 

where we have again put h\oo = 0.7. 

As a result, in the framework of the general STT, present cosmological observations, taking into account 
the Solar-system data, restrict the possible variation of G to values less then 10 _13 /yr. This estimate may 
be considerably tightened if the matter density parameter f2 m and the (negative) deceleration parameter q 
will be determined more precisely. 



G 



3 G-dot in (l+3+7V)-dimensional cosmology with multicomponent 
anisotropic fluid 

We consider here a (4 + 7V)-dimensional cosmology with an isotropic 3-space and an Einstein internal space 
[3"Tl 143j . The Einstein equations provide a relation between G/G and other cosmological parameters. 

3.1 The model 

Let us consider (4 + iV)-dimensional theory with the gravitational part of the action 

S 9 = ^j d i+N x^R , (3.1) 
where k 2 is the fundamental gravitational constant. Then the gravitational field equations are 

^ = ^ (T M_ (5 M_r_ )! (32) 

where Tp is a (4 + 7V)-dimensional energy-momentum tensor, T — Tjjf , and M, P — 0, . . . , N + 3. 
For the (4 + iV)-dimensional manifold we assume the structure 

M A+N = R, x Ml x K N (3.3) 

where R* is 1-dimensional time manifold, M% is a 3-dimensional space of constant curvature, = 
S 3 , R 3 , L 3 for k = +1,0, —1, respectively, and K N is a JV-dimensional Einstein manifold. 
The metric is taken in the form 

ds 2 = g M Ndx M dx N = -dt 2 + a 2 (t)g(f {x)dx i dx j + b 2 (t)g^ (y)dy m dy n , (3.4) 

where i, j, k = 1, 2, 3; m, n,p = 4, . . . , iV + 3; g\j , gmn , a(t) and b(t) are, respectively, the metrics and scale 
factors for M% and K N . 

For Tp^ we adopt the expression of the multicomponent (anisotropic) fluid form 

m 

(Tp 1 ) = ]T diag(-p«(t),p£ {t)5)Mt)5™). (3.5) 

a=l 

Under these assumptions the Einstein equations take the form 

o- atV 2 m 

i + T = $> (iV + ~ 3%Q ~ (3 - 6) 



Q = l 



a 2a 2 TVdb 2fc 



a a 2 o& a 2 AT + 2 

a— 1 

6 3d6 A • 2 



^-^[ P « + (7V-lK-iV^], (3.7) 



Here 



^n[5 W ]=A ff W, (3.9) 
m,n = 1, . . . ,N, where A is constant. The 4-dimensional density is 

p a,(4) W = /• ^^W&"( t ) p «( t ) = f^fy (3 . 10) 



where we have normalized the factor b(t) by putting 



d N yVgW = 1. (3.11) 
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On the other hand, to get the 4-dimensional gravity equations one should put 8nG(t)p a ^ (t) — n 2 p a (t). 
Consequently, the effective 4-dimensional gravitational "constant" G(t) is defined by 

8irG(t) = n 2 b- N (t), (3.12) 

whence its time variation is expressed as 



G/G = -Nb/b. (3.13) 



3.2 Cosmological parameters 



Some inferences concerning the observational cosmological parameters can be extracted just from the equa- 
tions without solving them [31] . Indeed, let us define the Hubble parameter H, the density parameters f2 a 
and the " deceleration" parameter q referring to a fixed instant to in the usual way 

H = a/a, n a = 8ttGp q ' (4) /3H 2 = n 2 p a /3H 2 , q = -ad/a 2 . (3.14) 

Besides, instead of G let us introduce the dimensionless parameter 

g = G/GH = -Nab/ab. (3.15) 

The present observational upper bound is g j 0.1, if we take in accord with |36[ 138] 

G/G < lCT 12 (y -1 ) (3.16) 

and H = (0.7 ± 0.1) x lO" 11 ^" 1 ) » 70 ± W(km/s.Mpc). 

3.3 The vacuum model with two Einstein spaces 

Here we consider the vacuum case when Tp 1 — 0. Let us suppose that to is an extremum point of the 
function b(t), i.e. b(t ) = 0. In this point we get G(t Q ) = 0. From (|3~6|l . (|3~7]) . ([3~8]) we get that for t = t 

3a Nb , 

— + — =0, 3.17 
a o 

a 2a 2 2k 

- = T-S> 3 - 18 

a a A a z 

Let us suppose that we "live" near the point to, then according to modern observations on acceleration 
of expansion of the Universe ( |44[ 145) ) we should put a(to) > and a(to) > 0. This implies k < due to 
(]3 . 18[) and b(to) < 0, A > due to p,17p and (|3 . 1 9[) . Thus, our 3-dimensional space should have negative 
curvature and the internal iV-dimensional space should have a positive curvature. 

From (|3. 1 7(1 - (|3 . 19|) we obtain using the definitions of cosmological parameters 

' 2/ ' T + M, (3.20) 



H 2 a 2 



W = 3|»l- ( 3 - 21 ) 



Here ao = a(to) and bo = b(to). 

Since we suppose that we "live" now near the point to, then we get 



b _ bo(t-to) 

b b {S ' 
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and due to (]3 . 13[) and (|3 . 1 7|) we find 

G 



- = -Nb/b « -N^(t - to) = 3 " o( ^ to) . 
G a 



The subscript "0" refers to to- Using (|3.14p we obtain in our approximation 

-3q H%(t-t ). 



G 
G 



(3.23) 



(3.24) 



Remind that go < and hence G/G > for t > to and G/G < for t < to; in our approximation G/G does 
not depend upon the dimension of internal space N = dim A'. 



3.3.1 Exact 1 + 3 + 6 solution 

Now we consider the exact solution from Ref. [IB] defined on the manifold 

M = R* x M (3) x M (6) , 

with the metric 

ds 2 = (A/a)"* [-2/f 2 (dr) 2 + \\ 3 \g%\x)dx l dxJ + f 2 \X e \g^(y)dy m dy« 
where (M^\g<-^) and (M {G \g^) are Einstein spaces: 

Ric[g«] = A lff «, 



i = 3, 6. Here we use the notations: A3 = 2/c and Ag = A. 
In f3~26ll 



A = W 2 



£3 



A = -3e 6 (t 2 + e 3 ) [1 + t (h(r, e 3 ) + d)] + e 3 e 6 > 0. 
Ci is constant and Ei = sgn(Aj), i = 3,6, and 



h{r,e 3 ) = - In 



T — 1 



r + 1 

/i(t, £3) = arctan(r), £3 = 1 



£3 = -1, 



(3.25) 



(3.26) 



(3.27) 



(3.28) 
(3.29) 



(3.30) 
(3.31) 



As was mentioned above we should restrict our consideration to the case when our 3-dimensional space 
has negative curvature and 6-dimensional internal space has positive curvature, i.e. 



£3 = -1, e 6 = 1. 

The analysis carried out in [IB] tells us that the scale factor of our 3-space 

a3 = a=(AAr 1/4 |A 3 | 1/2 



(3.32) 



(3.33) 



has the minimum at some point when the branch of solution with t € (t_,t + ) is considered. Here t_,t + 
are roots of the equation A( r ) = belonging to the interval (0,1). In this case the scale factor of our 
space &3(t) is monotonically decreasing in the interval (r_, t*) and monotonically increasing in the interval 
(t*,t+). 

The scale factor of internal 6-space 



rt 6 



= &=(AA)- 1/4 / 2 1/2 |A 6 | 1 



/2 



(3.34) 
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has maximum at some point To. It is monotonically increasing in the interval (t_,to) and monotonically 
decreasing in the interval (to, t+). For other branches of solution with either r S (r_ , t+oo) or r G (—00, r+), 
(|r_|, \t + \ > 1) we get monotonic behavior of scale factors 03(7") and ae(r). 
Let us consider the behavior of our solution in the synchronous time 

ds 2 = —dt 2 + al{t)gf{x)dx i dx^ + a 2 (t)g£J(y)dy m dy n , (3.35) 

where 



V2 j T dr'tf^r^fi 1 . (3.36) 



The function t s (r) is monotonically increasing from t s (r—) = to T = t s (r_|_). 

The scale factor of 3-space has minimum at to = t{ T *)- a 3(t) is monotonically decreasing from infinity to 
finite value in the interval (0, to) and monotonically increasing to infinity in the interval (to,T). 

The scale factor of 6-space has maximum at = <(t*). The function a§{t) is monotonically increasing 
from infinity to finite value in the interval (0,i*) and monotonically decreasing to infinity in the interval 
(t*,T). Only in the case when C± > we get that i* < to and hence in the time epoch near to we get an 
accelerating expansion of our 3-space. 

3.4 The model with two Ricci-flat spaces and two-component fluid 

Here we consider another example when two factor spaces are Ricci-flat. 
In this case, excluding b from (|3.6|) and (|3.8|) . we get 

N - 1 m 

-j^g'-g + q-Y,^^ ( 3 - 37 ) 

with 

A a = [2N + 1 + 3(1 - N)vt + 3iV^], (3.38) 

where 

vt=V%lp a , v a N =VNl9 a , P a >0- (3-39) 
When g is small we get from (|3.3T|) 

g^q-Y,A a n a . (3.40) 

Q = l 

Note that (|3.40p for N = 6, m = 1, v\ = v\ = (so that A 1 = 13/8) coincides with the corresponding 
relation of Wu and Wang [57] obtained for large times in case k = — 1 (see also [30] ). 

If fc = 0, then in addition to ([3T0]) . one can obtain a separate relation between g and Q , namely, 

7V — 1 m 

-^T5 2 -5 + i-E^ = ° ( 3 - 41 ) 

a=l 

(this follows from the Einstein equation R® — — k 2 Tq, which is a linear combination of (|3 . 6[) - (|3 . 8[) . 
3.4.1 Two-component example: dust + (N — l)-brane 

Consider now two component case: m = 2 03]. Let the first component (matter) be a dust, i.e. 

v\ =1/^ = 0, (3.42) 
and the second one (quintessence) be a (JV — l)-brane, i.e. 

v\ = 1, 4 = -1. (3.43) 
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We remind that, as it was mentioned in [48], the multidimensional cosmological model on the product 
manifold R x M\ x . . . x M n with fields of forms (for review see [14]) may be described in terms of multi- 
component "perfect" fluid [19] with the following equations of state for as component: pf = — p a if p-brane 
worldvolume contains Mi and pf — p a in the opposite case. Thus, the field of form matter leads us either 
to A-term, or to stiff matter equations of state in internal spaces. 

In this case we get from (|3.40[) for small g 

g^q-— — ^ +4— — -Q 2 , (3.44) 
y N + 2 N + 2 ' v ' 

and for k — and small g we obtain from (|3.4ip 

l-g^^ + fl 2 . (3.45) 
Now we illustrate these formulas by the example when N — 6 (K 6 may be a Calabi-Yau manifold) and 

- q = Q 1 = Q 2 = 0.5. (3.46) 

We get from ([3~44]) 

-0.06 (3.47) 

in agreement with (??). 

In this case the second fluid component corresponds to magnetic (Euclidean) A^55-brane (in D = 10 type 
I, Het or II A string models). Here we considered for simplicity the case of the constant dilaton field. 

This example tells us that for small enough temporal variation of G we may find the estimates on G-dot 
without consideration of exact solutions. But here we should select the solutions that give us the accelerated 
expansion of our world. We may use, for instance, the mechanism suggested above: but instead of curvatures 
in 2 factor-spaces we should consider two fluid components. 



4 Multidimensional cosmology with anisotropic fluid |49L 150] 
4.1 The model 

Now we consider a cosmological model describing the dynamics of n Ricci-flat spaces in the presence of 
1-component " perfect-fluid" matter [19] . Metric of the model 

n 

g= -e?q>[2~f(t)}dt(g) dt + ^2exp[2x l (t)]g\ (4.1) 

i=0 

is defined on the manifold 

M = Rx Mo x ... x M n , (4.2) 

where manifold Mi with a metric g l is a Ricci-flat space of dimension dj, i = 0, . . . , n; n > 2. The multidi- 
mensional Hilbert-Einstein equations have the following form: 

K-^R = K 2 T^f, (4.3) 

where k 2 is the gravitational constant, and the energy-momentum tensor is adopted as 

(T#) = diag(-p )Pl <J£\ . . . , Pn 6™ : ). (4.4) 

describing anisotropic fluid in general. 

We put pressures of the anisotropic perfect fluid in all spaces to be proportional to the density 

1 1 ■ 

Pi (t) = (l--^)p(t), (4.5) 
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where it, = const, i = 0, . . . ,n. Here we put p > 0. 

We impose also the following restriction on vector u = (m) E R n 

<u,w>*<0. (4.6) 

Here bilinear form < ., . >*: R n x R n — > R is defined by the relation 

<u,v >*= G^UiVj, (4.7) 

M, v e J? n+1 , where 

(5^ 1 

are components of matrix inverse to a matrix of the minisuperspace metric [511 123] 

Gij = d i 5 lj - didj. (4.9) 
In llUg) D = 1 + ^™ =0 is the total dimension of M P~2"j) . The restriction (l4"?rJ|) reads 

4.2 Solutions with power-law and exponential dependence of scale factors 

Let us consider two special families of solutions |19l 120] with the metric written in the synchronous time 
parametrization 



g = -dt s ®dt s + ^af{t s )g\ (4.11) 



i=0 

These solutions have either power-law or exponential dependence of scale factors (w.r.t. f s ). 

4.2.1 Solutions with power-law behavior 

The solutions with power-law behavior of scale factors take place for 

< u (A) -u,u 0. (4.12) 

The vector = 2di corresponds to the A-term fluid with p; = — p ( vacuum- like ) . In this case the 
solutions are defined by the metric (|4.1ip with the scale factors a, = a,(i s ) = Atf^ , and the density 

2 -2 < u,u >* 

Kp= <u^-u,u>itr (4 - 13) 

< u {Ji > — u, u >* 

with u l = G lJ Ui and Ai are positive constants, i = 0, . . . ,n. 

We will use the following explicit formulas for contravariant components 

1 " 

u' = G«ty = ^ + — ^tf, (4.15) 



eL 2 - £> 



and the scalar product 



< - W "«,«>.= - £ + ^ £ «, + ^ £ «0 2 - (4.16) 



i=0 i=0 ~ i=0 
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4.2.2 Solutions with exponential behavior 

The solutions with exponential behavior of scale factors take place for 

< u (A) - u,u >*= 0. (4.17) 

In this case, the solutions are determined by the metric (|4.11[) with the scale factors a,; = di(t s ) — 
A4 exp(iy l t s ), and the density p = const, 



S = £U i IJlfP, (4.18) 

where e = ±1, u l = G^Ui and Ai are positive constants, i — 0, . . . , n. Here p > for (u,u)# < and p < 
for (u, it)* > 0. 

Remark. The model under consideration was integrated in [19] for < u,u >*< 0. The solutions from 
[19] were generalized in [20] to the case when a massless minimally coupled scalar field was added. Families of 
exceptional solutions with power-law and exponential behavior of scale factors in terms of synchronous time 
variable were singled out in [20j and correspond to a constant value of the scalar field: ip = const. When the 
scalar field is omitted we are lead to solutions presented above (in [Tj|] these solutions were originally written 
in the harmonic time parametrization) . It may be verified that the exceptional solutions with power-law 
dependence of scale factors are also valid when the restriction (|4.6j) is omitted. Moreover, it may be shown 
that for < u,u >*= the power-law solutions are coinciding with our vacuum Kasner-like solution. In this 
case the matter source vanishes since: p = in (|4.13l) . 

4.3 Acceleration and variation of G 

The subspace (Mo,g°), g° to be flat, do = 3, describes our 3-dimensional space and (Mi,g l ) internal factor- 
spaces. We are interested in solutions with accelerated expansion of our space and small enough variations 
of G obeying experimental tests at the present moment [33] 



G 



(t s0 ) < 0.1, (4.19) 



GH 

Here we suppose that internal spaces are compact. Hence our 4-dimensions constant is (see [3Tj ) 

n 

G = const n« d >). (4.20) 

i=l 

4.3.1 Power-law expansion with acceleration 

For solutions with power-law expansion the accelerated expansion of our space takes place for 

v° > 1. (4.21) 
For D = 4 when internal spaces are absent we get 

u° = , (4.22) 
6 - u 

<u (A) -u,u>*= i(w -6)u t^O (4.23) 
6 

that implies uq ^ and uo 7^ 6 (here < u,u >*= — ^Uq < 0). The condition v° > 1 is equivalent to 
4 < uq < 6 or, equivalently, 

-P<P<-|. (4-24) 



13 



that agrees with the well-known result for D = 4. 
For the power law solutions we get 



G ELi TT a i'- 



.0 



=-!! = _. (4.25) 
and hence 

The constant parameter <5 describes variations of the gravitational constant and according to (|4.19|) 
|<S| <0.1. 

It follows from the definition of v l in (14. 14|) that 



^=\ Ud \ (4.27) 



or, in terms of covariant components (see (|4.15[) ) 

. (D-4)« -2E^ 1 « i 



i(5-£>) U0 +Er=i' 



(4.28) 



Thus, the relations (|4.16[) , (|4.21[) . (|4.28[) and the constraint (|4.12p define a set of parameters u t compatible 
with the acceleration and tests on G-dot. 

In what follows we will show that these relations do really define a non-empty set of parameters u» 
describing equations of state. 

The case of constant G. Consider now the most important case 6 — 0, i.e. when the variation of G is 
absent: G = 0. Indeed, there is a tendency of lowering the upper bound on G, e.g. according to arguments 
of [52] 6 < 1CP 4 . This constraint just follows from the identity 

G a 

that may take place in some multidimensional models (a is the fine structure constant). 

Isotropic case. First, we consider an isotropic case when pressures in all internal spaces are coinciding. 
This takes place when 

Ui = vdi, (4.30) 
i = 1, . . . , n. For pressures in internal spaces we get from (|4. 5|) 

Vl = (1 - v)p, (4.31) 

i = 1, . . . , n. In the isotropic case we get from (14. 10)) and (|4. 

<u,u >*= 2 ~D ^~3^~ + 2du ° v ~ 2dy2 }> ( 4 - 32 ) 

< u (A) -u,u >*= - -^ [2m + 2du + i(d - l)u% - 2du v + 2dv 2 }. (4.33) 
For S — we get in the isotropic case v = ^ or in terms of pressures 

Pi = |( 3 Po - 1, (4.34) 



Substituting into (TQ2")) and (|Q5|) we get 



< u, it >* = - ~Uo, (4.35) 
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< m (A) -u,u = ^u (u - 6). (4.36) 

Remarkably, we obtain the same relations as in D = 4 case (see Remark above). For our solution we 
should put Mo 7^ and u ^ 6. 

Using (|4.30p we get u° — —uo/6 and u 1 = for i > 0, hence z/j = for z = 1, . . . , n, i.e. all internal spaces 
are static. 

Metric (|4.11|) reads in this case 

n 

g = -dt s di s + 4jt»" ° 5 ° + £ A\g\ (4.37) 

8=1 

where are positive constants, and 

= (4.38) 
We see, that the power v° is the same as in D = 4 case. For the density we get from (|4.13[) 

* 2p= (^W (4 - 39) 



Thus, equations of state (|4. 5|) with relations (|4.30[) imposed, lead to the solution (|4.37[) - (|4.39|) with Ricci- 

p 

3" 



flat (e.g. flat) 3-metric and n static internal Ricci-flat spaces. For 4 < uq < 6, or — p < po < — § , we get an 



accelerated expansion of our 3-dimensional Ricci-flat space. 

Non-isotropic case. Let us consider the anisotropic (w.r.t. internal spaces) case with 5 = 0, or, 
equivalently (see (|4.28p ) , 



This implies 



(D-4)u = 2^Ui. (4.40) 

i=l 

< u (A) - u,u >*= -UoK - 6) - A, (4.41) 
6 

< u, u >*= --Uq + A, (4.42) 



G 



where 



2 

> 0. (4.43) 



no 1 / n \ 
i=l 1 \i=l / 

The inequality in (|4.43[) could be readily proved using the well-known Cauchy- Schwartz inequality: 

>(f>«) a . (4.44) 

i— 1 i— 1 i—l 

Indeed, substituting 6j = v5 an d Q = Ui/s/dl into (|4.44[) we get (|4.43|) . The equality in (|4.44[) takes 
place only when the vectors (bi) and (cj) are linearly dependent that for our choice reads: Uj/ \fd[ = v\fd[ 
where v is constant. Thus, A = only in the isotropic case (|4.30|) . In the non-isotropic case we get A > 0. 

In what follows we will use the relation 

< u (A) -u,u >*= -(n - u+ )(u - Uq), (4.45) 

where 

u± = 3 ± V9 + 6A (4.46) 
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are roots of quadratic polinomial (|4.41[) . obeying u < 0, Uq > 6 for A > 0. It follows from (I4.40[) that 
u = —uq/6 and hence 

"'- Mu^-eA - "»"*■ ( " 7) 

The function v°(uo) is monotonically increasing: i) from to +oo in the interval (-co, Uq); ii) from — oo 
to +oo in the interval (v,q,Uq); hi) from — oo to in the interval (uq,+oo). 

The accelerated expansion of our space takes place when v > 1, or, equivalently, when 

{A) u € (u^Uo), (4.48) 

(B) u €«„<). (4.49) 
ujfc = 2 ± V4 + 6A. (4.50) 



In terms of wo-parameter: 



p =w p, w = l-y, (4.51) 



these two branches read: 



U> /r :) ■ = v /1 + |a< w < l + ~Jl + ^A = w^ (4.52) 



(B) «,+ = -yi + ^A<«j <i-^yi + |A = «;+. (4.53) 

The first branch (A) describes a superstiff matter (wo > 1) with negative density. Indeed, the relation 
p < follows from (|4.13[) and < u,it >*> 0, see (|4.40l) . 

The second branch (B) corresponds to matter with a broken weak energy condition (since wq < — ^) and 
positive density (since < u, u >*< 0). This matter is a fantom one (i.e. w < —1) when A > 2. For A < 2 
the interval (wq,Wq^) contains both fantom points (wq < —1) and non-fantom ones (wq > —1). 

The case of varying G. Now we consider another important case <5 ^ 0, i.e. when the variation of G is 
non-zero: G ^ 0. We use the bound |<5| < 0.1, stating the smallness of 6. Using (|4.28p we get 

n 1 

^2ui = -dbu , (4.54) 

i=l 

where d — D — 4 and 

1 4. 1=4$ 

b = b(5)= 7 ™ ■ ( 4 - 55 ) 
For the scalar product we get from (|4.54[) 



1 2 



< u (A) -u,u >* = - Bu - A, (4.56) 
< u,u >*= -— ?ig + A, (4.57) 



where A was defined in (f4~4^|) (see (|4TTU|) and P~Td]) ). and 



A_A W >_,--g + iS'l, (4.58, 



12d(l-f) 2: 

1 - 



2 ' 

3 



S = B ( 6 ) = ( 4 - 59 ) 

1 2 

It should be noted that due to \S\ < 0.1 A is positive 4 > and close to 1: \A— 1| < il0~ 2 . 
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For contravariant component u° we get from (|4.15l) and (|4. 54|) 



u° = ~uo, (4.60) 
C = C(5) = '3B -2= T -^-g. (4.61) 



!-2 



It follows from (|4~56)l and g60]) that (see (|4~T4]0 

2Cu 



- Sit - A 
Here uq ^ u q , where 



^° = —a " r- ( 4 -62) 



= ^(5) = ^-(35 ± + 6AA) (4.63) 

are roots of quadratic polinomial (|4.56[) . 
In what follows we will use the identity 

n Aul - 4u Q - 6A , 

v - 1 = — r^ 2 2 r- 4.64) 

Au 2 - 6Bu - 6A y J 

Isotropic case. Let us consider an isotropic case (|4.30[) . In this case we obtain from (|4.54|) 

v = ^dbu . (4.65) 

or, in terms of pressures 

Pi = \[ib PQ + (2 - 3b)p],i = 1, . . . , n. (4.66) 



For scalar products we get 

< n. ii. > = — 

6 



<u,ii>«= — -ug, (4.67) 



< u (A) -u,u >* = ^-uoK - 6B). (4.68) 



For our solution we should put uq ^ and uq 6B/A. 
The metric (I4.11[) reads in our case 
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-dt s ® dt s + 4^"V + e g A 2 gij (4 69) 



1=1 



where A$ are positive constants, 

^ = (4 J0) 

Am - 65 v 7 
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d(l- |)(A«q-6B) 
i = 1, . . . , n. The last formula follows from (|4.14p and 

u< = — -5 

6d(l- §)' 

We see, that the power i/° is not coinciding for 5 ^ with that from D — 4 case. 



(4.71) 



(4.72) 
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For the density, since A > 0, we get from (|4. 13[) 



12A 



The condition of accelerated expansion of 3-dimensional space i/° > 1 reads as 



4 6B(S) 

< it < ' (4.74) 



A(S) ^ A(S) 
or, equivalcntly, in terms of wo-psLTamctcr (po = wop) (|4.51[) 

w ° (6) = 1 " ^iw <w ° <1 - = (475) 

For S > we get an isotropic contraction of total internal space Mi x . . . x M n . In this case tUg" (S) < — 1 
and hence phantom matter may occur with the equation of state close to the vacuum one, since 

w+(5) + l = ^1±J) , (4.76) 



For small 5 we have w£ (S) = 1 - f + 0(5 2 ). 

For 6 < we get an isotropic expansion of total internal space. In this case Wq (6) > — 1 and the phantom 
matter does not occur. In both cases u>o*(<0 < — ^ and Wq^(5) + 3 = 0(5 2 ). 



Non-isotropic case. Now we consider a non isotropic case A > when 5 7^ 0. 
Using relation (|4.64| we obtain 

(up - U^)(UQ - UQ t ) 

(uo - u£)(u - Uq) 

where = u^(S) were defined in (|4.43|) and 



»o» = " 



£ (<S) = 2 ± V4 + 6AA. (4.78) 

The accelerated expansion of 3-dimensional space takes place when v° > 1, i.e. 

(A) «oe(ii ,M,« (i)), (B) u G (ix+(<5),u+(<5)). 

In terms of u>o-parameter po = w^p, (wq — 1 — ^f-) these two branches read: 

(A) Wo(5)<w <w^(6), (B) w+(S)<w <w+JS), (4.79) 

where 

= <^) = 1-^. (4.80) 

For small 5 we have 

t4(J) = <(0) - 2(1 ± ^==) + 0(<5 2 ), (4.81) 

t4*(*)=^*(0)+O(* 2 ). (4.82) 
Thus, for small S the lower and upper bounds on wo have a small deviation from those obtained in the 
case 5 — 0. For small 5 the upper bounds shift only on 0(S 2 ) term while the lower bounds shift on 0(S) 
term. 

The first branch (A) describes a superstiff matter wq > 1, since Wq(S) > 1 due to Uq (5) < 0. It may be 
shown that the density is negative in this case since < u, u >*> 0. 

For the second branch (B) we get for the upper bound Wq^(S) < — 1/3 due to Uq^(S) > 4. For the lower 
bound we find that Wq (S) < — 1 only if 



A>6(A(5)-B(6)) = -^—^. (4.83) 

This is the condition on appearance of the phantom matter. For 5 > this inequality is valid, but for 
S < it is satisfied only for big enough A. 
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4.3.2 Exponential expansion with acceleration 

For solutions from 4.2.2, an accelerated expansion of our space takes place for v>° > 0. For D — 4, when 
internal spaces are absent, we get u° = —Uq/Q and (u,u)* = — Aug, (u^ — u, = ^(u — 6)u — 0, which 
implies uq — 6, or, equivalently, p = —p. 
We get 

V ° = ~ £ \/?' (4 - 84) 
which agrees with the well-known result for D = 4 de-Sitter solution with cosmological constant A = n 2 p > 0. 
The condition v° > is equivalent to e = — 1. 
For our exponential solutions we get 

% =-£V*. H=^ = V °, (4.85) 

and hence 

n 

G/{GH)= ^i/'dj =(5, i.e. (4.86) 

we get the same relation as in (|4.26p . 

The constant parameter 5 describing variation of the gravitational constant, obey the restriction \S\ < 0.1. 
Due to (|4.86[) we get the same relations (|4.27l) and (|4.28[) for S as in the powerdaw case. 

The case of constant G. Isotropic case. Consider the important case 5 = 0, i.e., when the variation 
of G is absent: G = 0. 

First, we consider the isotropic case when pressures coincide in all internal spaces, see (I4.30p . Here, we 
obtain the same relations as in D = 4 case. For our solution, we should put uq ^ and hence, due to (I4.17p . 
m = 6, i.e . po = -p. 

Using (|4.30[) . we get u° = — ito/6 = — 1 and u l — for i > 0, hence z/,- = for i = 1, . . . , n, i.e., all internal 
spaces are static. 

The metric (I4.11[) reads in this case as 

n 

g = -dt s ®dt s + Alex.p(2v°t s )g° + J^g\ (4.87) 

i=i 

where Ai are positive constants, and 

-° ^ -e\j^f- (4.88) 
For accelerated expansion we get e = — 1. We see that the power v° is the same as in D = 4 case. 

Anisotropic case. Consider now the anisotropic (w.r.t. internal spaces) case with 5 = 0, or, equivalently 
when (|4.40p is satisfied. It follows from (|4.40[) that u° = — uq/6 and hence 



"0 / 12K 2 P ; 

v =^y V / ^A> "o = . - (4-89) 



The accelerated expansion of our space takes place when v° > 0, or, equivalently, when either 

(A) u =Uq, e = l or (B) u — , e = -1. (4.90) 
In terms of the parameter wq, see ()4.51|) . these two branches read: 



(A) wo =w = Jl + fA, (4.91) 



(B) w = w+ =- x /l + |A. (4.92) 
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The first branch (A) describes the super-stiff matter (wq > 1) with negative density p < 0. 

The second branch (B) corresponds to matter with positive density (since (u, u)* < 0). This matter is 
the phantom one (i.e., wq < — 1) when A > 0. 



The case of varying G. Now, we consider the case S ^ 0, i.e., when G ^ 0. We take the observational 
bound | S | < 0.1. Using ([4728]) . we get relations ([454]) and ([455]) . It follows from (14756]) and ([460]) that 



n Clin / 12re 2 p , , „„. 

1/0 = " £ ^V^A- (4 ' 93) 



Here uq = Wq (^) are defined in (|4.63|) . 



Isotropic case. Let us consider the isotropic case (|4.30[) . We should put uq ^ and hence uq = 
6B/A > 0. The metric (|4~TT]) reads 

n 

g = ~dt s ® di s + exp(2z/ t s ).g° + exp(2^t s ) ^ A\g\ (4.94) 

i=l 

where are positive constants, 



Cu I2n 2 p 

u =-e—\\^r, and (4.95) 



i__ /___*_, ,-=!,...,„. (496) 



6d(l — <5/ 2 ) V Au o 

We see that the power v° does not coincide for 5 ^ with that in D = 4 case. 
The accelerated expansion condition for 3D space, v° > 0, reads as 



6B(S) 

no = ^f, e = -l (4.97) 

or, equivalently, in terms of wq (|4.51| (po = wop) 

-u,,s 2B((5) 

For 5 > 0, we get an isotropic contraction of the whole internal space Mi x . . . x M n . In this case 
Wq(5) < — 1 and hence phantom matter occurs with the equation of state close to the vacuum one since 

For small S we have Wq(S) = — 1 — 8/3 + 0(S 2 ). 

For 5 < we get an isotropic expansion of the whole internal space. Then, Wq(S) > — 1, and phantom 
matter does not occur. 



Anisotropic case. Consider the anisotropic case A > when 5^0. Here uq — Uq(S), see (|4.63|) . 
Accelerated expansion of 3-dimcnsional space takes place when v° > 0, or, equivalently, when cither 

(A) u = Uq(5), e = l or (B) u = u+(S), e = -I. (4.100) 

In terms of the parameter wq {po = wqp, wq = 1 — -y-) these two branches read: 

(A) w = wo (6), (B) w = wt(S), (4.101) 
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where Wq(S) = 1 — Uq(8)/3. 

For small 8 (see (|4.81[) the parameter Wq(8) has a small deviation from that obtained for (5 = 0. 

The branch (A) describes super-stiff matter wq > 1 since Wq (6) > 1 due to Uq(S) < 0. It may be shown 
that the density is negative in this case since (u, it)* > 0. 

For branch (B) we find that Wq(S) < — 1 only if (|4.83[) is satisfied. 

This is the condition on appearance of the phantom matter. For 8 > this inequality is valid, but for 
8 < it is satisfied only for a big enough value of anisotropy parameter A, see (14.831) . 



5 5-brane solution with acceleration and small variation of G 
5.1 The model 

In this section we deal with S*-brane solutions describing two electric branes and a set of I scalar fields 
The model is governed by the action 



S = J d D x^\g~\{R[g] - Kpg MN d M y a d N ^ 

E -^exp[2A a M](F a ) 2 }. (5.1) 

ra=l,2 a ' 

Here g — gMN{x)dx M ® dx N is a metric of pseudo-Euclidean signature (—,+,...,+), F a = dA a is a 
form of rank N a , (h a p) is non-degenerate symmetric matrix, <p — ((p a ) S R is a vector of / scalar fields, 
X a ((p) — X aa ip a , is a linear function. Here a = 1, 2 and a, ft = 1, . . . , I. In (|5.ip \g\ = | det(gMN)\- 
We consider as an example the manifold 

M = (0, +oo) x Mi x M 2 x M 3 x M 4 x M 5 . (5.2) 

where Mi are oriented Ricmannian Ricci-flat spaces of dimension di, i = 1, . . . , 5, and di = 1. 

Let two electric branes be defined by sets h = {1,2,3} and I2 = {1,2,4}. They intersect on Mi x Mi. 
The first brane covers also M3 and the second one covers M4. The first brane corresponds to the form F 1 , 
and the second one corresponds to the form F 2 . 

For world-volume dimensions of branes we get d(I s ) — N s — 1 = l+^+^+s, s = 1, 2 and d(Ii<~)l2) = l+c?2 
is the dimension of branes intersection. 

Consider now S'-brane solution governed by the function H = 1 + f p 2 , where p is a time variable, 
P = XQ 2 /8, 8 = 1,2. 

K = K S = d(I s )(l + ^L) + \ sa \ sP h a P, (5.3) 

s = 1, 2 is supposed to be non zero. Thus, K\ — K2 — K. Here (h a/3 ) = (hap)^ 1 . 
The branes intersection rule is following one 

d{h n h) = + A lQ A 2 ,^ - ijf. (5.4) 

This relation corresponds to Lie algebra A2 [54l [14] . Remind that if s = (t/ s ,C/ s ), s = 1,2, where 
"electric" U s vectors and scalar products were defined in [2T1H~5] (see also [11)1154] ). Relations K\ = K 2 and 
(|5.4p follow just from the formula (A as /) = (2(C/ S , U S ')/(U S ' ,U S ), where (A as /) is the Cartan matrix for A 2 
(with A 12 = A 21 = -1). 

We consider the following exact solution 

.9 = P 2A {-dp ®dp + H-^ip'g 1 + g 2 ) + PT 2 V + P" 25 .? 4 + g 5 } , (5.5) 

exp(^ tt ) =P SA ° +BA ?, (5.6) 
P 1 = -QH- 2 pdp A n A t 2 A 73, (5.7) 
P 2 = -QH- 2 pdp A n A T 2 A 7-4, (5.8) 
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where 



s=l,2 

B = 2R-\ (5.10) 

s = 1,2. Here g\ = dx®dx, t\ — dx and t; denotes a volume form on Mi. Remind that all Ricci-flat metrics 
<? 2 , . . . , g 5 have Euclidean signatures. 

This solution is a special case of a more general solution from [35] corresponding to Lie algebra A%, 

5.2 Solutions with acceleration 

Let us introduce a "synchronous" time variable r = r(p) by the following relation: 

P dp[H(p)] A (5.11) 

We put P < and hence K < that implies A < 0. Consider two intervals of the parameter A: 
(i) A < -1 and (ii) - 1 < A < 0. 

For the case (i) the function t = r(p) is monotonically increasing from to +oo, for p £ (0,pi), where 
pi = |P| -1 / 2 , while for the case (ii) it is monotonically increasing from to the finite value t\ — r(pi). 

Let space M5 be our 3-dimensional space with a scale factor 05 = H . 

For the branch (i) we get an asymptotical relation 05 ~ const t v , for t — > +00, where 

1 (5.12) 



-4+1 



and v > 1. For (m) we obtain as ~ const (ri — r) iy , for r — > n — 0, where < due to (|5 . 1 2|) . Thus, 
we get an asymptotical accelerated expansion of 3-dimensional factor space M5 in both cases i) and ii) and 
0,5 — > +00. 

Moreover, it may be readily verified that the accelerated expansion takes place for all r > 0, i.e. 
a'5 > 0, d'5 > 0. Here and in what follows we denote / = df /dr. 
Indeed, using the relation dr/dp — H A (see (JSTTTJ)) we get 

. dpda 5 2\A\\P\p 

a5 = ^ = ^r~ ' (5 - 13) 



and 



ft = £^ = S(l + |JV). (5-14) 



dr dp dp H 2+A 
that certainly implies inequalities for derivatives of scale- factor a§. 

Now we consider the variation of G. For our model the 4-dimensional gravitational coupling (in Jordan 
frame) is 

G = const ■ JI*_ (o t ~*) = H 2A P~\ (5.15) 

where 

ai =H A - 2B Pl a 2 =H A - 2B , a 3 = a± = H A - B (5.16) 

are scale factors of "internal" spaces Mi, . . . , M4, respectively. 

The function G(r) has minimum at the point tq corresponding to 

At this point the variation of G is zero. This follows from explicit relation for dimensionless variation of G 

1 - \P\p 2 

6 = G/(GH)^2 + ^ J U^, (5.18) 
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where H — ^ is the Hubble parameter. The function G(t) is monotonically decreasing from +00 to 
Go = G(t ) for r e (0,ro) and monotonically increasing from Go = G(r ) to +00 for r S (r ,fi). Here 
fi = +00 for the case i) and fi = n for the case ii). The scale factors 02 (r), ^(t), a 4 (r) are monotonically 
increasing from 1 to for r S (0, fi), since the powers A — B and A — 2i3 are positive and P < 0. The scale 
factor ai(r), is monotonically decreasing from zero to ai(r 2 ) for r £E (0,T2) and monotonically increasing 
from ai(T2) to zero for r £ (T2,fi), where T2 is the point of maximum. 

We should treat only solutions with accelerated expansion of our space and small enough variations of 
the gravitational constant obeying the present experimental constraint 

\8\ < 0.1. (5.19) 

Here like in case of the model with two curvatures [56 r is restricted by the interval containing tq. It follows 
from (|5.18[) that in the asymptotical regions S — ■> 2 that does not agree with experimental bounds (|5.19p . 
This restriction is satisfied for the interval containing the point To where 6 = 0. 

The calculation of G-dot in the linear approximation near To, gives the following approximate relation 
for dimensionless parameter of reciprocal variation of G [53j 

<5« (8 + 2|^|- 1 )Uo(t-t ), (5.20) 

where Hq — H(tq) (compare with the analogous relation in [56)). This relation gives approximate bounds 
on values of time variable r allowed by the restriction on G-dot. 

It should be stressed that the solution under consideration with P < 0, di = 1 and rfs = 3 takes place 
when the configuration of branes, the matrix {h a p) and dilatonic coupling vectors A a , obey the relations 
(|5.3|) . (|5.4|) with K < 0. This is not possible when (/ic/j) is positive definite, since in this case K > 0. In the 
next section we give an example of a setup obeying (|5.3p and (|5.4p . by introducing "phantom" fields. 



5.3 Example 

Let us consider the following particular example: iVi = N2, i.e. the ranks of forms are equal, and I = 2, 
(h a p) = —{5 a p), i.e. there are two "phantom" scalar fields. We also put = d±. 
Then relations (15.31) and (HL4|) read as 



and 



Xl = \l=d(I)(l + ^L)-K, (5.21) 

AiA 2 =d n + ffi^ + ^, (5.22) 

where d(I) = d(Ii) = d(I 2 ) = 1 + d 2 + d 3 , d n = d(I x n I 2 ) = 1 + d 2 and K < 0. Relations (j5T2T|) and (fST2"2"|) 
are compatible since it may be verified that they imply 

A|Aj e(-l,+i), (5-23) 



|Ai||A 2 

i.e. vectors Ai, A2 belonging to Euclidean space R 2 and obeying relations (|5.2ip . (|5.22p do exist. The left 
side of (|5.23p gives cos 9, where 8 is the angle between these two vectors. 
We get in this special case 

For K — > —00 the allowed time interval (t_,t+) of accelerated expansion obeying G-dot restriction (|5.19p 
vanishes, i.e. t + — r_ — > (sec (|5.20p ). 
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6 Conclusions 



In this paper we considered different cosmological models in diverse dimensions leading to relatively small 
time variation of the effective gravitational constant G. 

We estimated the possible variations of the gravitational constant G in the framework of a generalized 
(Bergmann-Wagoner-Nordtvedt) scalar-tensor theory of gravity on the basis of the field equations, without 
using their special solutions. Specific estimates were essentially related to the values of other cosmological 
parameters (Hubble and acceleration parameters, dark matter density etc.), but the values of G/G compatible 
with modern observations did not exceed 10 2 . 

We considered also the multidimensional cosmological model with an m-component anisotropic ("per- 
fect") fluid. The multidimensional Hilbert-Einstein equations led to relations between G and cosmological 
parameters. 

In case of two factor spaces with non-zero curvatures without matter, we have suggested a mechanism 
for predicting small G. When the 3-space has a negative curvature and the internal space has a positive 
curvature, we got at some time interval an accelerating expansion of our 3-dimensional space and a small 
value of G/G. We have shown that this result is compatible with the exact 1 + 3 + 6 solution from 46J. 
(Recall that only three exact solutions are known for a vacuum cosmological model with a product of two 
Einstein spaces, see [46]). 

We also presented another example where two factor spaces are Ricci-flat and for a two-component 
example (dust + 5-brane) we obtained a small enough variation of G. 

Besides, we considered multidimensional cosmological models describing the dynamics of n + 1 Ricci-flat 
factor spaces Mi in the presence of a one-component anisotropic (perfect) fluid with pressures in all spaces 
proportional to the density: pi — wtp, i — 0, n. Solutions with accelerated expansion of our 3-dimcnsional 
space Mo and small enough variation of the gravitational constant G were found. These solutions have 
either exponential or power-law behavior of scale factors w.r.t. synchronous time variable. In both cases 
they exist for two branches of parameter wq. The first branch describes superstiff matter with wq > 1 (and 
negative energy density) , while the second one may contain phantom matter with wq < — 1 (and positive 
energy density). Here, contrary to the two-curvature model, the experimental bounds on G are satisfied for 
all allowed values of the synchronous time variable. 

We considered an S'-brane solution with two non-composite electric branes and a set of / scalar fields as 
well. The solution, corresponding to Lie algebra A2, contains five factor spaces, and the fifth one, M5, is 
interpreted as our 3D space. As in the model with two non-zero curvatures, we found that there exists a 
time interval where accelerated expansion of our 3-dimensional space co-exists with a small enough value of 
G/G obeying the experimental bounds. Similar results for other rank 2 algebras were obtained in |57|. 

Thus, here we have shown that there exist different possible ways of explanation of relatively small 
time variation of the effective gravitational constant G compatible with modern cosmological data (e.g. 
acceleration) : we may consider either 4-dimensional scalar-tensor theories or multidimensional cosmological 
models with different matter sources. The experimental bounds on G may be satisfied either in some time 
interval or for all allowed values of the synchronous time variable (from (0, +00) for power-law case or from 
(—00, +00) for the exponential case. 

We considered recently [58] also the multidimensional gravity with a Lagrangian containing the Ricci 
tensor squared and the Kretschmann invariant. In a Kaluza-Klein approach with a single compact extra 
space of arbitrary dimension, with the aid of a slow-change approximation (as compared with the Planck 
scale), we built a class of spatially flat cosmological models in which both the observed scale factor a(r) and 
the extra-dimensional one, 6(r), grow exponentially at large times, but 6(r) grows slowly enough to admit 
variations of the effective gravitational constant G within observational limits. 
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